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This paper is concerned with exponentially ill-posed operator equations with additive 
impulsive noise on the right hand side, i.e. the noise is large on a small part of the 
domain and small or zero outside. It is well known that Tikhonov regularization with an 
data fidelity term outperforms Tikhonov regularization with an fidelity term in 
this case. This effect has recently been explained and quantified for the case of finitely 
smoothing operators. Here we extend this analysis to the case of infinitely smoothing 
forward operators under standard Sobolev smoothness assumptions on the solution, i.e. 
exponentially ill-posed inverse problems. It turns out that high order polynomial rates 
of convergence in the size of the support of large noise can be achieved rather than the 
poor logarithmic convergence rates typical for exponentially ill-posed problems. The 
main tools of our analysis are Banach spaces of analytic functions and interpolation- 
type inequalities for such spaces. We discuss two examples, the (periodic) backwards 
heat equation and an inverse problem in gradiometry. 
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1 Introduction 

In this work we analyze Tikhonov-type regularization for exponentially ill-posed problems where 
the data are corrupted by impulsive noise. We suppose that the measurements are described by 
functions ^ G (M) on a submanifold M C where g"^ denotes the exact data and 

^ the noise function. The noise ^ is called impulsive if |^| is large on a small part of M and small 
or even zero elsewhere. Impulsive noise naturally occurs in digital image acquisition due to faulty 
memory locations or any kind of physical measurements with malfunctioning receivers. Inverse 
problems with impulsive noise have been studied extensively in the literature (see e.g. [2H1] the 
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references therein), and several authors [T1f51f5|r7lllOII26] also analyzed Tikhonov-type regularization 
where reconstructions are defined as minimizers of a generalized Tikhonov functional 


fa e argmin 
f&D(F) 


1 


I(/)-<?' 


obs I 


^(/) 


( 1 ) 


Here df is a Banach space, F : D{F) C X ^ (M) the forward operator, TZ : X ^ (— 00 , 00 ] a 

convex penalty term stabilizing the reconstructions, r > 1, and a > 0 is a regularization parameter. 
Most common examples for IZ are powers of Banach space norms TZ{f) = ^ \\f — foWx with fo & X 
and g > 1 or total variation-type functionals. 

It is well known from several numerical experiments [5l ll41IT5] that the choice r = 1 in (CD leads 
to much better reconstructions than r = 2 , and several efficient algorithms have been developed 
to minimize CD with r = 1 (see [inilTHllinillllCZKlH] ) ■ The reconstruction improvements are even 
more striking if the forward problem arises from parameter identification problems in PDEs with 
smooth solutions [3]. 

The aforementioned analysis of CD always requires IICIIl'-(m) 
explain the remarkable difference between r = 1 and r = 2. Recently, a new model to describe the 
impulsiveness of ^ has been proposed by the last two authors in |13| . Following this approach, we 
will assume that 

3 P G 05 (M) : ll^ll< e, |P| < rj- (2) 

Here 05 (M) denotes the Borel tr-algebra of M and e, 77 > 0 are noise parameters. Under this 
model, the main result of [T3] can be described as follows: If F maps X Lipschitz continuously 
into a Sobolev space (M) with k > d/p and if the smoothness of the exact solution /I is 
described by a variational inequality with index function p (cf. (jH) for details; this is ensured 
under standard conditions as we will discuss after Assumption [T] and in Corollary 12.31) , then the 
Bregman distance between fl and fa is bounded by 


V 


{fcP 


< Cl— + C 2 --h C3{—ip)* (- 


( 3 ) 


with the Fenchel conjugate (—g?)*(s) := supj>Q(st + (p{t)) of —p (extended by 00 on (— 00 , 0)). If 
we set 77 = 0 and hence e = this result basically coincides with the estimates proven in 

the literature on regularization with Banach norms [niSlfSlTTlIlO) . Nevertheless, for mostly impulsive 
noise we expect e ~ 0 and very small 77 > 0. In this situation the remarkable improvement is given 
by the higher exponent 2 -I- l) — - of 77 , which is fully determined by the smoothing properties 

of F. ^ 

In this work we study CD in case that the range of F consists of analytic functions. This situation 
typically arises e.g. in parameter identification problems in PDEs with remote measurements. 
Thus we will obtain ([3|) for arbitrary k and p (probably influencing the constants), implying 
super-algebraic decay in 77 . Intuitively, as the smoothing property of the forward operator is of 
exponential type (if measured by the decay rate of the singular values), we may hope for an even 
better expression in 77 , namely an exponential one. Taking into account that in exponentially 
ill-posed cases the function tp is typically only logarithmic, this then implies logarithmic rates in 
e, but still polynomial rates in 77 under an optimal choice of a. 

This study is organized as follows: In the following § 2 we recall and adapt general convergence 
rate results for Tikhonov regularization CD under variational source conditions to our noise model 
in combination with a new interpolation-type inequality. The variational source condition will be 
validated under classical spectral source conditions. § 3 contains the construction and properties of 
Banach spaces (M) of analytic functions on three manifolds M (circles, intervals, and spheres) 
as well as certain interpolation-type inequalities for these spaces. The index A is always a weight 
function characterizing the growth of analytic extensions. If the forward operator additionally 
maps Lipschitz continuously into (M)) then the general convergence analysis from § 2 is ap¬ 
plicable. This mapping property is then verified for two practical examples in § 4, the (periodic) 
backwards heat equation and an inverse problem in gradiometry. We end this paper with some 
conclusions in § 5. 
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2 Generalized Tikhonov regularization 

This section is devoted to the analysis of Tikhonov regularization (ED- For the whole section let 
X he a Banach space, r > 1 and F : T>{F) G X ^ U’ (M) an operator. 

As discussed in |13] a minimizer /„ of ([T]) exists under reasonable assumptions on F and TZ. 
Furthermore, it can be proven by standard arguments that fa is stable w.r.t. 5 °'^® in a suitable 
sense and converges to as ||^||l'-(m) —>■ 0 if a = a (||^||l^(m)) is chosen such that ||CllLr(M[)/Q! —>■ 0 
and a —>■ 0. For details we refer to [131 Sec. 2.1] and the references therein. 

2.1 Assumptions and merit discussion 

In the following we state and discuss our assumptions used to prove rates of convergence. As 
usual an variational regularization theory will establish convergence estimates w.r.t. the Bregman 
distance V ^/q,, defined by 

V [fa, P) := (/) - (/i) - (r, / - p), 

where /* S dTl{p) is a subgradient of TZ at p. Note that 'D{f,P) = ||/~/^||^ if ^ (/) = 
11/ — /oil^ and A is a Hilbert space. 

Moreover, as in many other recent papers we will use an abstract smoothness assumption in the 
form of a variational inequality: 

Assumption 1 (Variational source condition). We assume that there exists /3 > 0 and a eoncave, 
inereasing function ip : [0, 00 ) —>■ [0, 00 ) with (/3(0) = 0 such that 

f3V{f,P)<n{f)-nP) + ‘p{mf)-F{p)\\lr^M^) forallfGV{F). (4) 

If F : A — 7 > (M) is a bounded linear operator between Hilbert spaces, TZ{f) = \\f — /o||^, 

fo G X and r = 2, then it has been shown that o is both sufficient and necessary for 

convergence rates of order ^||CIIl 2 (m)) ■ Note that this furthermore shows that a spectral source 
condition must imply a variational source condition with a function p yielding the same rates of 
convergence in this situation. Another proof of this implication has been given by Flemming [^ 
for general p. 

We are mainly interested in exponentially ill-posed problems, where the most relevant case for 
the function p in Assumption [1] is p(t) = Cpp{t) with 

'Pp(f) = t<exp(-l). (5) 

This function can be defined on all of K by concave extension, but for our asymptotic studies only 
the behavior close to t = 0 is of relevance. As shown by Flemming |^, if A is a Hilbert space, r = 2 
and F is bounded and linear, the spectral source condition p = pp (F* F) w implies a variational 
source condition ([1]) with p = (3'p2p and /3' > 0 depending on ||r(;||, ||F||, and p. 

For nonlinear operators F, the condition Q can be seen as a combination of source and non¬ 
linearity condition (see e.g. m Props. 3.35 and 3.38]). For inverse medium scattering problems 
it was recently shown with the help of complex geometrical optics solutions (cf. [ 12 ]) that a varia¬ 
tional source condition with p of the form Q holds true if the solution belongs to a Sobolev ball, 
and p is an explicit function of the Sobolev index. 

Both techniques for verifying Assumption |T] mentioned above yield ((3|) with r = 2. To obtain 
L^-variational source conditions from L^-variational source conditions, we will use the following 
smoothing properties of F, which also play a central role in the rest of our analysis: 
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Assumption 2 (smoothing properties). Let X be a Banaeh spaee, F : ^{F) C A —>■ L^(M) and 
IZ . X ^ (— 00 , oo] a eonvex, lower semi-continuous penalty functional. Assume that there are 
(5o > 0, <7 > 1 and an inereasing funetion 7 : [0 , Jq] —>■ [0 , 00 ) with lim 5 _>o 7 (< 5 ) = 0 such that 

\\F{f) - 5^|Uoo(m) < y{S)V (/, + i||F (/) - (g) 

for all 0 < S < Sq and all f G X. 

Remark 2.1. Suppose there exists a Banach space Z C L°“(M) fl L^(M) such that 

(A) F maps Lipschitz continuously into Z, i.e. 

\\FU)-F{f^)\\z<CV{f,f^Y (7) 

with C > 0, q > 1 and 

(B) there exists i5o > 0 and an increasing function 7 : [0, Jq) —?> [0, 00 ) with yiS) = 0 such 

that for all 0 < S < 60 and all g G Z 

I|5||l“(m) < 7 ( 1^)115112 + ^II5||li(m)- (8) 


Inserting © in ([S]) with g = F {f) — gl proves that Assumption[^ holds true with y{d) in ([0]) 
replaced by Cy{6). 

An inequality of the form ([5]) is shown in m Lem. 3.2] in the case of Sobolev spaces on 
a bounded Lipschitz domain M C Z = (M) with k > -. In this case we obtained 

k 1 

7 (( 5 ) = ci6'^~p where ci > 0 is a constant depending only on M,/c and p, but not on 5 S (0,(5o). 
As (HI) holds for any 0 < <5 < (5o, it can (in this special case of 7 ) be reformulated as 


llslli 


< inf 
0<5<5o 


fe _ 
Ci0'^~ 




:||5ll 


LH 




with u := k/d — 1/p. If — ciitdo I|5 |Iw'^-p(m)> follows from calculating the infimum 

using differentiation w.r.t. 5. However, since IH^’^ (M) is continuously embedded into L^(M), ([S]) 
holds true for all g, with a larger constant C if ciuSq is smaller than the embedding constant. 
Consequently, is of the form of a classical interpolation inequality and ([5]) can be interpreted 
as an interpolation inequality as well. 

The more general formulations ([5]) and m allow for a wider class of spaces Z, in particular 
spaces of analytic functions, which are more suitable for exponentially smoothing operators F. 
As in the asymptotic behavior of 7 (5) as ^ 0 will depend on the smoothing properties of 

F, more precisely the more smoothing F the faster 7 (^) \, 0 as 5 0. In § 3 we will introduce 

spaces Z which guarantee (HJ with exponentially decaying 7 . 

Now we are in position to prove the following proposition which validates Assumption [1] under 
a spectral source condition. As we are mainly interested in exponentially ill-posed problems, we 
restrict ourselves to logarithmic source conditions, which are appropriate in this case. 

Proposition 2.2. Let X be a Hilbert space, M C and choose 72. (/) = ||/||^, r = 1. Suppose 
that F satisfies Assumptions^ with q = 2 and (HJ holds true. 

1. For all a > 0, 0 <r] < 5 q/ 2, and e > 0 the estimator /„ from ® is bounded by 


A' 


< \\f 


lx 


^ 27?7 ( 277 ) 


( 10 ) 
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2. Suppose satisfies the variational source condition (Assumption[Jj) with r = 2, a is chosen 
such that 

2 r]y {2r]) 


+ < Ce (11) 

a a 

for all e,rj with some constant Ce, and F satisfies the Lipschitz condition ©• Then there 
exists Cv > 0 independent of a, e and rj such that fl satisfies the variational source condition 
o with r = 1 for all f = fa and ip replaced by ip{Cv-)- 


Proof. Part 1: Using the minimizing property 
1 
a 


F{fa)-g°'^^ 
and the noise model ([5]) we find the estimate 


+ 'F{fa) < —\\£,\\l^{U) PF-iP) 
Li(M) a ' ' 


F{fa)- 9 


obs 


Li(M\P) 


a 


Using ([B]) with 5 = 2ri and V (^fa, = \\fa — P\\ < ||/a|| + ||/^|| (since X is a Hilbert space) 

it follows 


fo 


a' 


<£ + a\\P\\t+ F{fa)-g^ 


Li(P) 


( 12 ) 


(13) 


F{fa) - 9^ 
which then implies 

Fifa)-g^ 


< - 

Li(P) 2 


Fifa) - P 


Li(M) 


+ 777 


(2^)( 


fo 






< 


Li(P) 




As 


F{fa) - 9' 


obs 


LI 


F{fa) - 9^ 


LI 


(14) 


< e, we find from inserting m into (fTBl) that 


< ?£ + ||;t||^ + ( £ +||/.|| ). 

a " a V A* " "^/ 


277 (27?) 


A* a 


fo 

Now we add ^ 'fpP on both sides, complete the squares and obtain 


a' 


777 (277) 


<-+(ll/'l 

a 


777 (277) 


< 


- + l|/'l 

a " ' 


777 (277) 


a 


Taking square roots yields the desired estimate. 

Part 2: Note that \\F (f) - < H-F'(/) - 5licx,(M) ||^ (/) “ 51 li(m) ' Furthermore ® 

together with o yields 

\\f if) — gPi^ac(yi^ < 1 iP \\f — P\\x P \\f if) — 

< + ll/-/'IU 

for arbitrary 6. Together with the first part this shows that \\F {f)—g^\\p(^M) <CP\F{f)-pU i(M) 
with Cv := (7 ((5) + C/6)(2\\p\\x P Ce), which yields the assertion. □ 


As mentioned after spectral logarithmic source conditions imply variational logarithmic 
source conditions. As (p 2 p{C„x) = ip 2 p{x) P o(l) as a: \ 0, we obtain: 

Corollary 2.3. If X is a Hilbert space, F : X ^ L^(M) is linear, the spectral source condition 
P = (pp{F*F)w is satisfied, and (HU holds true. Then the variational source condition m holds 
true for f = fa with r = 1, and ip = fi'p 2 p where /3' > 0 depends only on Ill’ll, p, 5 q, and 

Ce. 
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2.2 Rates of convergence 

Now we will present error estimates for the minimizers fa in O- Our analysis uses the observation 
by Grasmair |10] that the approximation error can be bounded by the Fenchel conjugate function 
of —If. Let us introduce a function ip by 

ipia) := {-(fi)* f--) = sup ((p{t) - , a > 0. (15) 

V C(J r>0^ OiJ 


Note from this definition that ip is concave and non-decreasing. 

For ipp defined in ([5]) one obtains ipp{a) = ^In (1 -|- o(l)) as a —>■ 0. 

We are now in a position to formulate and prove the main result of this section, which closely 
follows [131 Thm. 3.4]: 

Theorem 2.4. Let X be a Banach space, F : I?(F) C X ^ (M) be an operator and IZ : X ^ 

(— 00 , 00 ] a convex, lower semi-continuous functional. Let fa be a solution of o with r = 1 and 
suppose that the noise function f fulfills ([2]). If Assumptions[J\and\^ hold true with g > 1 (actually 
o is only needed for f = fa), then for any b > 1 there exists rjo > 0 such that 


PV {Jaj^) 


a 


-777 


2b 


id 


Tl|-^(/a) ~ + 2 


pii ilB!’)) 

(,5a)9'-i 


6-1 

Q' 

-h 2aq''tfp{2ba) 


q'ip{ba), 


(16) 

(17) 


for all 0 < r] < r]Q, e > 0 and a > 0. Here and below q' denotes the conjugate index of q, i.e. 



Proof. By assumption there is a measurable P such that ||CIIli(m\p) < e and |P| < 77. Thus we 
have for any g G L^(M) that 

\\g — 5^IIli(m) + II?IIli(m) = Wg — 5^IIli(m\p) + \\g ~ + II^IUhM) 

< (II 5 ~ + IICIlLbM\P)) + II 5 ~ 5^IIl1(P) + IICIIlI(M) 

^\\g — + 2 ||^||li(m\p) + 2 ||g — g^\L^(p) 

<ll5 ~ + 2e -I- 2|P|||g — g^||Loo(M[). (18) 


For g = F{ja) we will now apply the interpolation inequality ([5]) with 5 := Given 6 > 1 we 

set 770 := and obtain for any 0 < 77 < 770 that 


1 


^ \\Fifa) - g'|| l 1(M) + IICIUbM) 

< ll^(/a) ~ g°'^'*llLi(M) + 2 ( e -b 777 


26 

6-1 


V^'D(fa,pyy (19) 
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Now Assumption [T] implies 

{LP) 


(gji 


< TZ{P) - nP) + (\\F{P) - F{p)UiiM)) 

g>{\\F{l)-FiP)\\ 


GU 1 1 - 

< —IICIIlI(M)- \\F{fa)-g' 

a a 

US 2 

< - 


iLi(M) 


Li(M) 


< 


+ ^{\\F{f^)-FiP)\\mM)) 


( 20 ) 


The second term on the right-hand side can be handled by Young’s inequality xy < x'^ /q' y'^/q 

with X = which implies that (ITBl) holds true. 

{L,p)>o, we derive from (HOD and (fTBl) in a similar fashion that 


Since V 


in 2 

< - 


26 

e + yy { 7 —-y V 


2 ba 


6-1 


( 7 a,/^) 


26a 


\\FiU)-p\\LHM) 


< 


2 

< -e 


+ y(ii2"(7.)-a/*)iiL.(M)) 

+ P (Pt") ^ (Ff')' 

+ 77^(s’"(s^'')) + V (/../')+v-(2H 

26 


GU 2 q' 


< 


jl-q' 


a 


-hi 

a 


6-1 


qtl){2ba), 


where we used that '0 is monotonically increasing. This immediately proves (dll). □ 

Remark 2.5. The optimal a can he obtained by balancing out the terms, cf. rm Thm. 2.3, 3.]. 
For simplicity we will do this only in the specific examples m § 

Note that the phenomenon of exact penalization (cf. will occur if ip (x) = c-x with some c > 0, 
as then tp (a) =0 if a < 1/c and ip (a) = oo otherwise. Consequently, in this situation the optimal 
parameter is a = 1/c, and for noise-free data this already yields an exact reconstruction. 

Remark 2 . 6 . All choices of a for which the right-hand side of (1161) tends to 0 (e.g. a-priori 
choices which approximately minimize this right-hand side) satisfy the condition ( 1111 ) . 


3 Spaces of analytic functions and interpolation 

inequalities 

Let A (M) denote the set of all analytic functions on M. In this section we will introduce spaces 
A''‘(M) C A (M) of analytic functions on the manifolds M = ]R/27rZ, (—1,1), and SS'^ and establish 
interpolation inequalities for spaces L°°(M) C L^(M) C A^(M). 


7 









Convergence Rates for Exponentially Ill-Posed Inverse Problems with Impulsive Noise 

General strategy 

We will introduce classes of analytic functions (M) which are indexed by an increasing function 
A : [0, oo) —>■ K U { 00 }. A will be a bound on the growth rate of analytic extensions of functions 
g G Al^(M). Moreover, for each choice of M we will introduce a sequence of finite dimensional 
subspaces Um C Al^(M), to G N with corresponding orthogonal projections Pm ■ L‘^{M) —>• Um 
and prove exponential bounds of the approximation error of the form 

Wg - Pm9\\L<^(U) < 9A(TO)exp(-A*(TO))||g||_4A(M) (21) 

for TO > Too and g G Al^(M) with the Fenchel-conjugate A*(s) := sup 2 ,>o[rs — A(r)] of the function 
A (extended by oo on (—oo,0)) and a polynomial g\. As a second ingredient we will need that 
the projection operators Pm can be extended to Pm '■ A^(M) —?> L°°{M.) and satisfy the operator 
norm bound 

ll-Pmff||L~(M) < Km||g||i,i(M) for all g G L^(M). (22) 

This is the case if Pm has an integral representation {Pmg){x) = Km(x, y)g{y) dy and || Arm||L“(Mx 

Km- The kernel Km is given by Km{x, y) = 9j{^)9jiy) fo^ orthonormal basis { 31 ,..., gNm} 

Um- With the help of (1^ and (1^ we obtain 

||5||l”(M) < Wg - Pmg\\L°°{U) + \\Pmg\\L°°{U) 

< gA(TO)exp(-A*(TO))||g|Ux(M) + KmllylUqM) 

for TO > Too and g G Al'^(M). We then choose m{S) (decreasing in d > 0) such that Km(S) < 

This yields the interpolation inequality 

I|5||l“(m) < 7(^)IIsIUmm) + ^IIsIIlTm) with 
7 ( 5 ) := qx{fn{6)) exp{-X*{m{S))) 

for all g G Al^(M) and (5 < do if 'm(So) = too (see dS]))- 

3.1 Case M = M/27rZ 

Note that any 27r-periodic real-analytic function g : K./27rZ —>■ R can be extended to a holomorphic 
function g on some strip S'b = {- 2 GC| |3(z)|<i?} with B > 0. We may measure the smooth¬ 
ness of g in terms of B and the growth of | 3 (a;-|-iy)| as |y| ^ B. This leads to the following 
definition (cf. [T71 Sec. 11] for the special case that A is an indicator function): 

Definition 3.1. Let A : [0,oo) —>• R U { 00 } he an increasing weight function with positive B\ := 
sup{r G [0, 00 ) I A (r) < 00 }. With the holomorphic extension g of g on B\ introduced above we 
define the space (R/27rZ) by 

A^ (R/27rZ) := jy G C'(R/27rZ) | g exists on Sb^, \\9\\a='W2-k7.) < 00 } > 

llfflU^(R/ 27 rZ) := sup [exp(-A(|9(z)|)) | 3 (z)|]. 
zCSb^ 

Theorem 3.2. Let A : [0, 00 ) —>■ R U { 00 } be non-decreasing with B\ > 0. 

1. The space A^ (R/27rZ) equipped with |M|_ 4 X(R/ 27 rZ) a Banach space. 

2. The Fourier coefficients g{n) of g G A^ (R/27rZ) satisfy 


\g (n)| < exp (-A* (jnj)) ||3lUx(R/2,rZ) > 


n Gh. 


(24) 
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3. If Pm denotes the L'^-orthogonal projection onto span {exp (in-) | |n| < m}, then there exists 
mo G N and a constant ca > 0 such that 


IK-f--Pm) < 2 cAexp(-A* (m)) \\g\\A>-{^/2^’L) 


(25) 


for all g G (R/27rZ) and m > mo- 

4 . There exists Jq > 0 such that the interpolation inequality 
and 

y{S) = CAexp (-X* 


holds true with M = R/27rZ 


for all g G Al^(R/27rZ) and all 0 < 6 < Sq- Here [ccj := supjn Gl n < x} for a; G R. 

Proof. Part 1: Obviously, (R/27rZ) is a normed space. Since any Cauchy sequence in this 
spaces converges uniformly on compact subsets of B\, the pointwise limit is again holomorphic, so 
it belongs to A^ (R/27rZ). This follows from the Cauchy integral formula and proves completeness. 
Part 2: Let g G A^ (R/27rZ) and denote by g the holomorphic extension oi g to Sbx- Since g and 
g(2Tr + •) coincide on R by the periodicity of g, it follows from the identity principle that g and 
g{2TT + •) also coincide on i-S- 9 is 27r—periodic. 

For 0 < R < B\ consider the annulus Kb := {w S C | exp(—i?) < |ry| < exp(i?)}. The map 
9 : Sb ^ Kb,0{z) := exp (iz) is holomorphic for any 0 < R < B\. The restriction of 9 to the 
rectangle Db := (z G S'/? | 5R (z) G (—tt, tt]} is bijective and 9 (—tt + iy) = 0 (tt + iy) = — exp(—y). 
Hence h:= go 9~^ ■ Kbx — !> C is continuous on the line segment {— exp(—y) | — R < y < R} and 
therefore holomorphic in the annulus Kbx - If follows that for all w = exp (ix) G SS^,x G (— 7 r, 7 r] 
we have 

h{w) = go 9~^{w) = g{x) = ^ g{n) exp {jnx) = ^ g{n)w^. (26) 

Being holomorphic in the annulus Kb^ i h has a Laurent series expansion which by the identity 
principle is uniquely determined by the series expansion (Uni) on SS'^, i.e. h{w) = 
for all u> G For any s G R with exp {—Bx) < s < exp (Ha) we have 


i)l = 


1 

27r* 


/ Hw) 


dw 


n+1 


< max \h(w)\ 

s S 


1 


But max|^,|^s \h{w)\ = maxQ(2)^in(i/jj) |5('Z)|) and hence with r := In (1/s) this implies 
\g{n)\ < max \g{z)\ •exp(-r |n|) < ||5lU>(B/27rZ) ’ exp(- [r |n| - A(r)]) 

for all 0 < r < B\. Optimizing in r proves (1241) . 

Part 3: The orthogonal projection is given by Pmg{x) = Z|n|<m s(^) ^ ^ 

Using ([M]) we obtain 


II 5 - -Pm5llL“(-77,7r) = sup 

aiG (—tTjTt) 


g{n) exp (i 

|n| >m 


mx) 


< 2||5lUA(K/2,rZ) (")))■ 


Recall that A* is always convex. As A(r) = 00 for r < 0, A*(s) = sup^>Q[sr—A(r)] is non-decreasing, 
and since Ha > 0 it is easy to see that there exists mo G N with A* (mo) > 0 and a := A* (mo +1) — 
A*(mo) > 0. Thus the convexity of A* implies a < A*(m+1) —A*(m) for all m > mp. Consequently, 


for any n > m> mo, we have A*(n) — A*(m) = + 1) “ ^*(/)) > (n — m)a. Therefore, 


j=m 


Y exp(-[A*(n) - A*(m)]) < Y exp(-a(n - m)) = . \ 

whenever m > mo. This shows (USD with ca := 2 /(exp(a) — 1 ). 

Part 4-' We follow our general strategy with Km{x,y) := ^ (Dirichlet kernel). 


‘ — 


2m+l 

277 ’ 


and m((5) = [f — ■ 


□ 
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3.2 Case M = (—1,1) 

The main idea is to extend a real-analytic function g on (—1,1) to a holomorphic function g on 
an ellipse with foci {—1,1} of the form 


£r = I X -\-iyGC 


x,y G 


cosh (r) sinh (r) 


2 — 


< 1 


r > 0. 


Definition 3.3. Let A : [0,oo) —>■ M U {oo} be an inereasing weight function with positive B\ := 
sup{rG[0, oo) |A(r)<oo}. We define the space {—1,1) by 

-4^ (-1,1) := G1) \ g exists on £b^, I|5IUa(_i_i) < ooj , 


llfflU^(-i,i) := sup 

0<t<Bx 


exp(-A(t)) sup \g{z)\ 

z^dSt 


Theorem 3.4. Let A : [0, oo) — >■ M U {oo} be non-decreasing with B\ > 0. 

1. The space (—1,1) equipped with the norm is a Banach space. 

2. The coefficients an{g) of g G A^ (~1) 1) with respect to the Chebychev polynomials satisfy 

|an( 5 )| < 2exp(-A* (|n|)) , n e N. (27) 

3. If Pm denotes the L'^-orthogonal projection onto the space of polynomials of degree < m, 
then there exists mg G N and a constant ca > 0 such that 


II {I -Pm) g \\< CA exp (-A* (to)) (28) 


for all g G A^ (~1) 1) and to > toq. 


4 . There exists do > 0 such that the interpolation inequality (1231) holds true with M = (—1,1) 
and 

7(d) = CAexp ^-A* { -1)1 (29) 

for all g G A^{—1, 1) and all d < do- 


Proof. Part 1: This can be shown as in the proof of Theorem 13.21 

Part 2: The following approximation argument is taken from Kress [T71 Thm. 11.7]. The function 
9 : Kbx £bxi ^{'ca) + on an annulus Kbx as defined in the proof of Theorem 13.21 is 

surjective and holomorphic (but not injective). 

Now let g G .4''(—1,1) with holomorphic extension g : £bx C. Define h : —>■ C by 

h{w) := 2g o 9{w). The holomorphic function h in the annulus Kbx can be expanded into a 
Laurent series h{w) = ^ ^Bx with coefficients 


— 


5 o h" + - 


drc 


J\w\—e-x.-p{r) \2 \ ^JJ ^ 

Substitution w = — shows that a_„ = a„. Therefore 


„rt+l ' 


For jwl = 1 we write w = exp (it) and obtain 


—B\ < r < B\, n G Z. 


h{w) = aoOn (w'^ -I , w G Kbx ■ 

nGN ^ ^2 


+ ^ ) = cos(nt) = T„(cos(t)) = T„ ( ( ic -f 
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where T„ is the n—th Chebychev polynomial. For z G £bx w G Kb^ with z = 9 (w) to 

obtain 

g{z) =go 9{w) = ^h{w) = Y + X! 

nGN 

The formula for a„ yields for any n G N and 0 < r < B\ that 

2 

- exp(rn) ledS “ 2 exp (A(r) - nr) (30) 

where we used the definition of ||•||_ 4 X(_l Optimizing in r yields (071) . 

Part 3: As P^g {x)=ao{g)+ Y,n=i°-n{ 9 )Tn{x) and ||T„||ioo(_i^i) = 1, we obtain 

OO OO 

\\g-Pmg\\L’^(-i,i)= XI ^ 2 X exp(-A* (n)) || 5 ||_ 4 X (_1 p . (31) 

71=771+1 71=771+1 

The sum may be bounded as in the proof of Theorem 13.21 

Part 4- We again follow our general strategy. A complete orthonormal system of Um = 
span{a:°, • • • , x™} C L^(—1,1) is given by {\/j + l/2pj '■ j = 0,..., m} with the Legendre poly¬ 
nomials pj. As ||Pj||loo(-i,i) = 1, the supremum of the kernel K{x,y) = + ^/‘^)Pji^)Pjiy) 

is bounded by Km = J2jLo(J + 1/2) = (m -|- 1)^/2. Hence, we choose m(S} := [+2+] — 1. □ 

3.3 Case M = SS^ 

Let SS^ := {x G : |x |2 = 1} denote the unit sphere. Recall that if Pm denotes the space of 
polynomials in x = (xi,X 2 ,X 3 ) G of degree < m, then the space of spherical harmonics Hm is 
defined by 


Rm '■= {P|SS 2 1 p G Pm harmonic and p is homogeneous of degree to} . 

Here harmonic means Ap = 0 and homogeneous of degree to that p(rx) = r"^p(x) for all r > 0 
and X G We have the following decompositions as orthogonal direct sums with respect to 
(•,-)l 2 (SS 2 ) (see e.g. [21 Cpt. 4]): 


L^(ss^) = ^ni, rm\ss^=^ni 

z=o z=o 


(32) 


One has dun'Hm = 2to-|- 1 and dimPm\ss^ = X]/lo(2^ + l) = + The orthogonal projections 

of L‘^{SS‘^) onto Rm and Pm\ss^ will be denoted by Qm and Pm, respectively. Spherical harmonics 
are closely related to Legendre polynomials Pm- Choose any orthonormal basis (Lj)^™}*'^ of Rm 
with respect to {■,-)l^(ss^)- The addition formula of the spherical harmonics (see e.g. [191 Thm. 
2]) states that ^^^pm{{x,y)) = yjix)Yj{y), i.e. 

{Qm9){x) = — - / Pm{{x,y))g{y)dy, x G SS'^ (33) 

47r 7552 

Let us introduce the averaging operator M : C'(5'5'^) —>■ C'(S'S'^ x [—1,1]) by 


{Mg){x,t) 


hv^ss'^\{v,x)=t}9i.y)‘^y^ i.i)> 

g{±x), t = ±l. 


Note that {Mg){x,t) is the average of / over a circle of radius ^/1 — P around x and that Mg is 
in fact continuous if / is continuous. 
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Definition 3.5. Let X : [0,c») —>• MU {c»} be an increasing weight function. We define the space 
A^iSS"^) by 

A\SS^) := {g e C{SS^) \ {Mg){x, •) S ^^(-1,1) for all x e 55^, ||g|U.( 552 ) < oo} , 
ll5lU^(SS2) := sup \\Mg{x,-)\\AH-i,i)- 

xGSS2 

Theorem 3.6. Let X : [0, cxd) —>■ R U {c»} be non-decreasing with B\ > 0. 

1. eguipped with the norm || • ||^a(ss 2 ) is a Banach space. 

2. For all g € A^{SS'^) we have 

-|— \ 

IIQmg|lL~(ss 2 ) < —— exp(-A*(m - 1))||5|U>(SS2). (34) 

3. There exist constants c\, c^a > 0 and mo € {0,1,...} such that 

Wg - f’m5llL“(SS2) < (CA + mdx) exp(-A*(m))|| 5 |UA(ss 2 ) (35) 

for all g G A^{SS^) and m > mo. 

4 . There exists 5o > 0 such that the interpolation ineguality (IMl) holds true with M = SS^ and 



for all g G A^(SS^) and all 0 < 6 < Sq. 

Proof. Parti: Obviously24^(5'S'^) is a normed space. To show completeness, let {gn) be a Cauchy- 
sequence in A^ISS"^). Then [MgA) converges to some G G C'(5'5'^x [—1,1]) with sup 3 ,gg 52 11 ^( 0 ;, OIU^c-i,!) < 
00 and {gn) converges uniformly to G(-, 1). Since M is continuous, it follows that G = MG{-, 1). 

Part 2: Due to ^ and the formula f{y) dy = | f(y) dydt we have 


{Qmg){x) = 


2 m + 1 


2 m + 1 


Pm{t){Mg){x,t) dt 


1-1 


Pm{cos{s)){Mg){x, cos(s)) sin(s) ds. 


With the mapping d{w) := ^^{w + ^) from the proof of Theorem 13.41 and the substitution w = 
exp (is) we obtain 


2m + 1 f 

{Q 7 ng){x) = --- K / Pm {0{w)) {Mg){x, 9{w)) dw. 

4 J{wec I |™|= 1 } 


We may deform the contour of integration {wj |w| = 1} to any contour {w| |r(;| = exp(r)} with 
|r| < B\. To estimate Qmg, we have to estimate the growth of the Legendre polynomials on £r. 
We use the identity Pm{z) = y Jq (z + \/z^ — 1 cos (p) d(p (see [ISJ Cpt. 4]) to find that 


Pm 



1 

TT , 
2 
TT 


0 L 


1 + cos ip 1 1 — cos ip 

W ---h 


W 


2 

-\ m 


dip 


w cos^ Ip -sim Ip 

w 


dip. 
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Therefore, \pm{z)\ = sup|^|^,, \pm {w + ^))| < exp(mr). It follows that 

\pm{0{w))\ \{Mg){x,e{w))\\dw\ 


\{Qmg)ix)\ < ^ sup 


0<r<B\ J |if;|=exp(r) 

p777 _|_ 1 

< —-— sup [exp (r(m - 1) - A(r)))] ||5lU^(ss=) 

4 r>0 


2 m + 1 


exp(-A*(TO - 1))||5 |Ua(ss2). 


Part 3: As shown in the proof of Theorem 13.21 there exists toq G {0,1,...} and a > 0 such that 
A*(mo) > 0 and A*(m + 1) — A*(m) > a for all m > mo- Using the identity 
for |r| < 1 and its derivative X)m=i = (1 ~ we obtain 

OO 

Wg - Pg\\L°°(ss^) < X! \\Qmg\\L°°{ss^) 


^ E 


2j + l 


exp(-A*(j - l))||g|U>(ss 2 ) 




< exp(-A*(m)) 


2j + l 


exp(-A*(j - 1) + A*(m))||g|Ux(ss 2 ) 




< exp(-A*(m))|| 5 |UA(s 52 ) ^ ^ ) exp(-a0 


1=0 


= exp(-A-(m))||9|U,„,., + 

This shows (ESI). 

PQiT^t J^. -By 0iHci dSSl) the kernel of the L^-orthogonal projection onto Pm\ss^ is given by 
Km{x,y) = ^ + As \\pi\\l°°(-i,i) = 1, the supremum norm of is bounded 

by Km = ^ + 1) = + 1)^- Therefore, we choose m{5) = [^dTr/AJ — 1. □ 

4 Examples 

In this section we show how our general techniques can be applied to practical examples. As we 
have seen in § 2, the key ingredients are the variational source condition o and the smoothness 
assumption dSD- As the variational source condition can be validated using Proposition 12.21 we 
will focus on verifying (|6|). This is done exemplary for the forward operators connected to the 
periodic backwards heat equation and an inverse problem in satellite gradiometry on SS^. Similar 
techniques should apply to many other exponentially ill-posed problems, e.g. inverse scattering 
problems. 


4.1 Backwards heat equation 

Let / G L^(—TT, tt). Let u : (— tt, tt) x [0, oo) —>■ M be a solution of the periodic heat equation 
^{x, t) = ^{x, t) if (x, t) G (-TT, tt) X (0, oo) 

u{a;,0) = f{x) if a; G (—7r,7r) (initial condition) (36) 

u{—TT,t) = u{Tr,t) if t G (0, oo) (boundary condition) 

It describes heat propagation on a circle parameterized by x as angular variable. The initial 
boundary value problem (1361) has a unique solution given by 


i{x, t) = exp f{n) exp (ina:), (a:, t) G [— tt, tt] x [0, oo) 


(37) 
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where f{n) are the Fourier coefficients of /. 

Fix a time i > 0. In the following we will study the forward operator F associated to the 
backwards heat equation, which consists in determining the initial heat distribution / from the 
measurement of the heat distribution u{x, t) at time t. This problem is known to be exponentially 
ill-posed, which can readily be seen from the decay of the singular values exp (—n^t) in ill). 

Let us therefore define T : L^(—7r,7r) —>■ L^(—tt, tt) by 

exp f{n) exp (inx ), x G [—tt, tt] (38) 

n^Z 


and set g := u{-,i). 

Lemma 4.1. Let T as in (1551) . set X := (— tt, tt) and TZ (/) = ||/||^ 2 (_^ 

1. Then T : X ^ (]R/27rZ) is bounded with 


= rG[0,oo). 


Furthermore Assumption\^holds true with q = 2 and 

jlS) < - 77 p=r-exp I — 

^ - exp(3i)-l ^ I 


TT 1 

S~2 


t . 


2. If /t G HP{ —TT, tt), then Assumption \7\ holds true with any bounded D{F), M = (—7r,7r), 
r = 1, some 13,13' >0 and 

T {t) = /3'ipp (r), T > 0. 


Proof. Part 1: First we estimate \\T\\x^A^(m/ 2 TrZ)- For a; G M, t > 0 we will use the identity 
: ^ exp (—(n -|- x)^t) < 

nGZ ^ n^Z 


—j=z ^ exp (—(n -|- x)^t) < —= ^ exp {—n'^t) < v^max < 1 


1 

Vt 


(39) 


which follows from Poisson’s summation formula (cf. [551 Cpt. 5, Thm. 3.1]). 

Applying Cauchy-Schwarz, supQ(j,)^j_^ | exp(mz)| < exp(|n|r) and (I59ll with x = we find 


sup |r/(.)|<^ exp(—n^t) exp(|n|r) 


i(z)=± 


n^Z 

1 




\X 


exp(A(r))||/|U ^exp (jnj - t'j 
nez ^ ^ 


The sum on the right hand side can be bounded as follows: 


nez 


Eexp(-(i«i-^-) i) <E 


nez ^ 


(- 3 ) 


exp I — (— 77 =) ^ + oxp I — 




< 2 exp(—< 2'/^max |l, t . 


nez 


This implies ||T||A 7 -».al^(R/ 27 rZ) ^ 2max{l,t The assertions will now follow from Remark 

nn The conjugate function of A is A*(s) = s'^t. By Theorem 13.21 we find that ([5]) holds true for 
all 0 < d < <5o: do > 0 sufficiently small with 


7(d) 


exp (3f) — 1 


( 

TT 

1 

lexpl 

_d 

2_ 


14 
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since toq = 1 and a = A* (mo + 1) — A*(mo) = 3t in the proof of Theorem 13.21 This yields the 
claim. 

Part 2: One readily calculates that 

(pp{T*T){w){x) - exp (inx), 


i.e. ||(/5p(T*r)M;||^2p(_^_^) ~ Thus/■!■ G HP{-tt,tt) is equivalent to/'!' = v?p /2 {T*T)w 

with w G L^(—7r,7r). Now Proposition 12.21 with M = (—7r,7r) yields the claim. □ 


Theorem 4.2. Let t > 0 and T be as in (IM)l and suppose p G with some p > 0. 

Suppose furthermore that a = a{rj, e) is chosen such that dill) holds true in the limit r]^e,a\i 0. 

Then there exists rjo > 0 and a constant C = C (t,T,fl^ such that for any noise function 
^ G L^(—TT, tt) fulfilling @ with 0 < 7? < ryoj e > 0, the following estimates are valid for Tikhonov 
regularization (H)) with TZ (/) = ||/||^ 2 (_^ r = !.• 

||7a - PWh <c(j^ + ^exp + (-ln(a))”^^ (1 + o(l)), (40) 

\\T{fa) - 7||li < C ^e+ ^exp + a (-In (a))"^^ (1 + o(l)), (41) 

as ri,e,a\ 0. The parameter a = d(e, if) can be chosen such that 

ll/a - PWh = O (max|(-ln(e))”^,ry2p|) , 

l|T'(/a) = O ^max|e,exp 

as e,p\ 0. 

Proof. As (fTTl) holds true, we know from Proposition 12.21 that (g]) can be applied with / = /„. 
From Lemma l4.II we know that we may apply Theorem 12.41 with (p (r) = /3Vp P) 


l{5) < 


8max{l,t-i/2| 

- hri —r^exp 

exp (3t) — 1 


^ ^ max{l,f 
~ exp {'it) — 1 


exp 


The function if in m is given by 



(-J/) (1 + 0(1)), 


5\0. 


if (a) = P' ^In ^ (1 + 0 (1)) < (/3')^ P (- In («)) ^ (1 + o (1)). 

Inserting these results into (HU) and (113 with q = q' = 2 and 6 = 2 yields (go]) and dnj. 

An optimal a can be chosen now depending on which of the first two terms in d40ll dominates. 
If the first term is larger than the second, we choose ai to balance the first and the last term in 
dMl), which yields 

- = (-In (e))-^ (1 + 0 ( 1 )), ai,e^0. 

OC\ 

On the other hand, if the second term is larger than the first, we choose a 2 such that the second 
and third term equal, which gives 

{-\n{a 2 ))~^ = t~P (l + o(l)), a 2 ,r]^ 0 . 

Taking the maximum of both cases yields the proposed bounds. As mentioned in Remark 1 2.6 1 the 
condition (ITT]) is satisfied. □ 
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4.2 An inverse problem in satellite gradiometry 

Let us assume, the earth is described by the unit ball B := {x £ : |a;| < 1}. In geophysics 

the value of the gravitational potential u on the surface of the earth is of interest as it contains 
information about the interior of the earth. If the value / of u on dB is given, u is the solution to 
the exterior boundary value problem 

{ Au = 0 in \ B, 

u = f ondB = SS^, 

|u(a;)| = 0{\x\~^) as \x\ —>■ oo. 


Using satellites it is possible to measure the second derivative of u in radial direction r = \x\ at 
some distance R> 1 from the earth, i.e. the rate of change of the gravitational force: 

uqq 2 

g = ttt on. nob 
ar‘^ 

The inverse problem of gradiometry then consists in estimating / given g (see [Hi Sec. 8.2] and 
references therein). Representing u by the Poisson formula for the exterior of -R, it can be seen 
(see e.g. i) that the forward operator is given by 




IS'S'2 


/SS2 


dR^ 


(m + l)(m 

m=0 


R 


■ 2 ) 


^m+3 


_i l-R 
\R~'^x - 

{Qmf){x) 


X G SS^ 


(42) 


with the orthogonal projection Qm onto the spherical harmonics of order m introduced in (llldll . 
This representation of T shows again that the problem to recover / from g is exponentially ill- 
posed. 

Lemma 4.3. Let T as in (11^ with R > 1, X = L'^{SS^) and TZ (/) = 11 / 11 ^ 2 ( 552 ). 

1. T is a bounded mapping from L^(S'S'^) to A^{SS'^) for 


A(r) 


—41n(i? — exp(r)), 
oo 


0 < r < ln(i?), 
else. 


2. Assumption\^holds true with 

^{ 5 ) = C5-5/2^-\/4V5-4 (43) 

with C independent of R and 6. 

3. If p G (SS^), then Assumption[l\ holds true with any bounded D {F), M = 55^, r = 1, 
some a, /?' > 0 and 


ip (r) = aV2p (t) , T > 0. 

Proof. Part 1: First note that for g^ G Rm and u G {0,1,...,} we have 

Pnit) f 


nl fl 

Pn{t)(Mgm){xN) dt = 


/-I 


1 27r'\/l J{x,y)—t 

Pui{x,y))gm{y) dy 


9m{y) dy dt 


/SS2 

47r 


“ 2m -b 1 ^^^ 3 m)ix) 
dlT 

= dn.m'Z ' T9m{x). 
2m -b 1 
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As + 1 } is an orthonormal basis of L^(—1,1), we find that {Mgm){x,t) = ^ 4 ^*^ 9m(x)■ 

Together with we obtain 


^ Pm{t) {m-I l){m 2 ) 

iMTf){x,t) = 2^ -- {Qmf){x). 


m —0 


Using the bound \pmiz)\ < exp(mr) for 2 ; £ d£r from the proof of Theorem Id.Bl and \\Qmf\\L°°{SS^) < 
{2m + l)||/||i 2 ( 5 S 2 ) we obtain 

Il'T'/IUuss^) < sup sup sup exp{-X{r))\{MTf){x,z)\ 

r>0 zCdSr xeSS^ 

^ /n Pm{z)im + l){m + 2 ) 

(m +l)(m + 2 )( 2 m+l)e’' 


- NZr> ."'AUf." 

<CR\\f\\L^. 


The last inequality follows by evaluating derivatives of the geometric series, and C is a generic 
constant independent of R. 

Part 2: The supremum of r !->■ sr — A(r) is attained if e’' = sR/ (4 + s), so 


exp(-A*(s)) 


i±£Vri±£^ 

4R J \ sR J 


< Cs'^R-^-^. 


Now (l4^ follows from Theorem 13.61 

Part 3: The proof is similar to the heat equation case using m Prop. 16]. There it has been 
shown that 

/t = (pp (t*T) w, w G ^ fl e HP ( 552 ) . 

Thus the spectral source condition is satisfied by assumption, and it follows from Corollary 12.31 
that dH) holds true with r = 1, and (p = I3'p2p for any bounded D{T). □ 

Theorem 4.4. Let R > 1,T as in and p G HP (S'S'^) with somep > 0. Suppose furthermore 
that a = a {rj, e) is chosen such that (11111 holds true in the limit ry, e, a 0. 

Then there exists 770 > 0 and a constant C = C (i?, T, /I) such that for any noise function 
f, G L^{SS'^) fulfilling ([2]) with 0 < V Si Voj e > 0, the following estimates are valid for Tikhonov 
regularization ([ 1 ]) with Rif) = II/IIl 2 ( 552 ), r = 1 .- 

ll/a ~/^IIl 2 (ss 2 ) ^ C'^—+ ^^7? ^R Ina) ^^(1 + 0 ( 1 )), 

r(/a) - T{P)\\lHSSZ <c(^e+ + a{- Ina)""^) (1 + o(l)), 

as 77, e, a \ 0. The parameter a = d{e, 77 ) can be chosen such that 

ll/a - PWhiSS^) = ^ (max I (-In 77 ^ 1 ) 

Wnh) - T{p)\\miss^^ = O (max{e,77^i?-V^}) 

as e, 77 \ 0 . 

Proof. As (imi holds true, we know from Proposition 12.21 that (|3]) can be applied with f = fa- 
From Lemma 14.31 we know that we may apply Theorem 12.41 with p (r) = P'p 2 p (ct) and 7 as in 
dMl). This yields the claimed error estimates. 

To obtain the asserted convergence rates, we may again distinguish if the e term or the 77 term 
dominates. Then we balance the dominating term with the pure a term to obtain a choice of a 
satisfying (dH). The maximum of both cases yields the claim. □ 
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5 Conclusions 

We have extended a recent approach for inverse problems with impulsive noise to the case of 
infinitely smoothing forward operators using standard Sobolev smoothness assumptions for the 
solution. Remarkably, one obtains high order polynomial rates of convergence in the size 77 of the 
corrupted domain, even though the underlying problem is exponentially ill-posed in the classical 
sense. We examined two exponentially ill-posed problems arising in PDEs and showed that our 
analysis can be applied. 

Our study gives rise to several further interesting questions: One concerns the extension of our 
analysis to more general domains in A possible strategy could be the use of local averages 
as introduced here for the case of SS"^. However, boundaries would cause technical difficulties, 
and the choice of approximating subspaces is not obvious. So far we have only studied a priori 
parameter choice rules, i.e. the solution smoothness characterized by Lp has to be known for 
choosing a. The popular discrepancy principle does not apply in our context (as only ||C||j;,i(m\p) 
is small, but P is unknown), and thus other a posteriori parameter choice rules would be of interest 
for practical applications. Other open questions include the optimality of the error bounds and an 
analysis of iterative regularization methods for nonlinear inverse problems with impulsive noise. 


Acknowledgement 

Financial support by the German Research Foundation DFG through subprojects A04 and G09 
of CRG 755 is gratefully acknowledged. 


References 

[ 1 ] M. Burger and S. Osher. Convergence rates of convex variational regularization. Inverse 
Probl, 20(5):1411-1422, 2004. 

[2] R. H. Chan, C. Hu, and M. Nikolova. An iterative procedure for removing random-valued 
impulse noise. IEEE Signal Proe. Let., ll(12):921-924, 2004. 

[3] C. Clason and B. Jin. A semismooth Newton method for nonlinear parameter identification 
problems with impulsive noise. SIAM J. Imaging Sei., 5(2):505-536, 2012. 

[4] C. Clason, B. Jin, and K. Kunisch. A duality-based splitting method for C-TN image restora¬ 
tion with automatic regularization parameter choice. SIAM J. Sei. Compub, 32(3):1484-1505, 
2010 . 

[5] C. Clason, B. Jin, and K. Kunisch. A semismooth Newton method for data fitting with 
automatic choice of regularization parameters and noise calibration. SIAM J. Imaging ScL, 
3(2):199-231, 2010. 

[6] J. Flemming. Generalized Tikhonov regularization and modern eonvergence rate theory in 
Banaeh spaces. Shaker Verlag, Aachen, 2012. 

[7] J. Flemming and B. Hofmann. A new approach to source conditions in regularization with 
general residual term. Numer. Eunet. Anal. Optimiz., 31:254-284, 2010. 

[8] J. Flemming, B. Hofmann, and P. Mathe. Sharp converse results for the regularization error 
using distance functions. Inverse Probl, 27(2):025006, 2011. 

[9] W. Freeden, F. Schneider, and M. Schreiner. Gradiometry—an inverse problem in modern 
satellite geodesy. In Inverse problems in geophysical applications (Yosemite, CA, 1995), pages 
179-239. SIAM, Philadelphia, PA, 1997. 


18 



Convergence Rates for Exponentially Ill-Posed Inverse Problems with Impulsive Noise 


[10] M. Grasmair. Generalized Bregman distances and convergence rates for non-convex regular¬ 
ization methods. Inverse Probl., 26(11):115014, 2010. 

[11] T. Hohage. Regularization of exponentially ill-posed problems. Numer. Funct. Anal. Optimiz., 
21:439-464, 2000. 

[12] T. Hohage and F. Weidling. Verification of a variational source condition for acoustic inverse 
medium scattering problems. Inverse Probl., 31:075006, 2015. 

[13] T. Hohage and F. Werner. Gonvergence rates for inverse problems with impulsive noise. SIAM 
J. Numer. Anal., 52(3):1203-1221, 2014. 

[14] B. Jin. A variational Bayesian method to inverse problems with impulsive noise. J. Comput. 
Phys., 231(2):423-435, 2012. 

[15] Q. Jin and L. Stals. Nonstationary iterated Tikhonov regularization for ill-posed problems in 
Banach spaces. Inverse Probl., 28(10):104011, 2012. 

[16] T. Karkkainen, K. Kunisch, and K. Majava. Denoising of smooth images using L^-fitting. 
Computing, 74(4):353-376, 2005. 

[17] R. Kress. Linear integral equations. Springer-Verlag, New York, 1999. 

[18] Y.-R. Li, L. Shen, D.-Q. Dai, and B. W. Suter. Framelet algorithms for de-blurring images 
corrupted by impulse plus Gaussian noise. IEEE Trans. Image Process., 20(7):1822-1837, 
2011 . 

[19] G. Muller. Spherical harmonics. Springer, 1966. 

[20] M. Nikolova. Minimizers of cost-functions involving nonsmooth data-fidelity terms. Applica¬ 
tion to the processing of outliers. SIAM J. Numer. Anal., 40(3):965-994 (electronic), 2002. 

[21] M. Nikolova. A variational approach to remove outliers and impulse noise. J. Math. Imaging 
Vision, 20(l-2):99-120, 2004. Special issue on mathematics and image analysis. 

[22] O. Scherzer, M. Grasmair, H. Grossauer, M. Haltmeier, and F. Lenzen. Variational Methods 
in Imaging. Applied Mathematical Sciences. Springer, 2008. 

[23] E. M. Stein and R. Shakarchi. Fourier analysis: An Introduction. Princeton University Press, 
2003. 

[24] E. M. Stein and G. L. Weiss. Introduction to Fourier analysis on Euclidean spaces, volume 1. 
Princeton university press, 1971. 

[25] G. Szego. Orthogonal polynomials. AMS, 1959. 

[26] E. Werner and T. Hohage. Gonvergence rates in expectation for Tikhonov-type regularization 
of Inverse Problems with Poisson data. Inverse Probl., 28(10):104004, 2012. 

[27] J. Yang, Y. Zhang, and W. Yin. An efficient TVLl algorithm for deblurring multichannel 
images corrupted by impulsive noise. SIAM J. Sci. Comput., 31(4):2842-2865, 2009. 

[28] W. Yin, D. Goldfarb, and S. Osher. The total variation regularized model for multiscale 
decomposition. Multiscale Model. SimuL, 6(1):190-211 (electronic), 2007. 


19 



